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Recent progress in the design and fabrication of artificial two-dimensional (2D) materials paves 
the way for the experimental realization of electron systems moving on plane fractals. In this 
work, we present the results of computer simulations for the conductance and optical absorption 
spectrum of a 2D electron gas roaming on a Sierpinski carpet, i.e. a plane fractal with Hausdorff 
dimension intermediate between one and two. We find that the conductance is sensitive to the 
spatial location of the leads and that it displays fractal fluctuations whose dimension is compatible 
with the Hausdorff dimension of the sample. Very interestingly, electrons in this fractal display a 
broadband optical absorption spectrum, which possesses sharp “molecular” peaks at low photon 
energies. 


Introduction. — A variety of exerimental protocols that 
can be used to create artificial two-dimensional (2D) lat¬ 
tices for electrons, atoms, and photons are nowadays 
available. Examples include schemes for creating artifi¬ 
cial honeycomb lattices [1], where a wealth of interesting 
phenomena have been observed such as Mott-Hubbard 
split bands [2], massless Dirac fermion behavior mod¬ 
ified by pseudo-electric and pseudo-magnetic fields [3], 
and photonic Eloquet topological insulating states [4]. 
In the case of ultracold atomic gases loaded in honey¬ 
comb optical lattices, recent progress has even led to the 
experimental realization [5] of the Haldane model [6]. 

In the solid state, a combination of e-beam nanolithog¬ 
raphy, etching, and metallic gate deposition [7-11] can 
in principle yield high quality two-dimensional (2D) pat¬ 
terns with arbitrary shape in semiconductor heterostruc¬ 
tures (such as GaAs/AlGaAs) hosting ultra-high mobil¬ 
ity 2D electron gases (EGs). Ultimately, these proce¬ 
dures yield an external potential landscape with the de¬ 
sired geometry that acts as a lattice of potential wells to 
trap electrons. The spatial resolution of these techniques 
can reach values of a few tens of nanometer or even be¬ 
low. Eurther improvements in spatial resolution can be 
obtained by bottom-up nanofabrication methods such as 
nanocrystal self-assembly [12]. These approaches allow 
to independently control the electron density and inter¬ 
site distances and to tune the interplay between on-site 
and nearest-neighbor repulsive interactions and single¬ 
particle hopping, opening the way to the observation of 
collective phenomena and quantum phase transitions in 
such artificial solid-state systems [2]. Synthethic solid- 
state quantum materials can also been created by uti¬ 
lizing scanning probe methods [3]. Here, suitably cho¬ 
sen molecules can be positioned with atomic precision 
on top of a substrate (such as Cu) with the aid of the 
tip of a scanning tunneling microscope (STM). Electrons 
confined in the substrate surface states and subject to 
the potential created by the deposited molecules can be 
probed via STM measurements [3]. 

These experimental achievements motivate the theo- 



FIG. I. (Color online) Panel a) The Sierpinski carpet after m 
iterations. The initial red unit is a n x n square lattice (n = 2 
in this work). This unit is replicated A/" = 8 times to obtain 
the green structure, which is £ = 3 times longer in each of 
the two spatial directions and contains an inner hole with no 
lattice sites. This operation is repeated to obtain the blue 
structure. For a number of iterations m ^ I we obtain an 
approximation of the Sierpinski carpet, a plane fractal with 
Hausdorff dimension du = log^A/" (— 1.89 in this figure). 
Panel b) Energy dependence of the density of states n(E) per 
lattice site for a Sierpinski carpet with m = 3. Panel c) and 
d) Two examples of Sierpinski carpets with leads (orange) 
attached in different positions. 
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retical investigation of complex 2D structures, with the 
aim of discovering novel transport and optical features 
which could enable or improve technological applications. 
In this work we present a theoretical study of the trans¬ 
port and optical properties of a 2DEG in a Sierpinski car¬ 
pet (SC), which is a self-similar 2D structure [13] shown 
in Fig. 1. The self-similarity of the SC is mathematically 
quantified by the fact that its Hausdorff dimension [13] 
dn (i.e. a generalization of the topological dimension) is 
between one (a line) and two (a plane), which makes the 
SC a fractal [14]. 

While brownian motion and the heat diffusion equa¬ 
tion on fractal geometries have been extensively studied 
in the literature [15, 16], the transport and optical prop¬ 
erties of electrons roaming on such complicated geomet¬ 
rical structures have comparatively received less atten¬ 
tion. More precisely, some analytical [17, 18] and nu¬ 
merical [19-22] studies of the conductance of electrons in 
Sierpinski fractals have appeared in the literature. We 
are instead unware of studies of the optical properties 
of electron systems in fractals, with the exception of an 
experimental work [23] on a one-dimensional potential 
realized by growing a fractal sequence of quantum wells 
in (Al,Ca)As heterostructures. 

In this work we show that the conductance of a 2DEC 
on a SC can reach the maximum value allowed by the 
number of open channels in the leads, depending on the 
lead positions and their widths, and it displays fractal 
fluctuations [24, 25] as a function of energy, in the absence 
of a magnetic field. We show that extended states, which 
are responsible for large conductance values, are quite 
robust to elastic disorder. Finally, we calculate the real 
part of the optical conductivity, which turns out to be 
much larger than e^/h and highly structured up to a 
photon energy of the order of t/2, where t is the hopping 
amplitude. 

We hasten to stress that the problem at hand is very 
different from that of a quantum particle displaying a 
self-similar spectrum. Such problems are very well stud¬ 
ied in physics, a paradigmatic example being that of the 
Hofstadter butterfly spectrum [26] displayed by an elec¬ 
tron moving in 2D under the combined effect of a periodic 
potential and a perpendicular magnetic field. Finally, 
we are not interested in the distribution of eigenvalues 
and nature of the corresponding eigenstates of electrons 
in plane fractals, which have been studied in great de¬ 
tail [27-31]. Rather, our aim is to unveil fundamental dc 
and ac transport characteristics, which can be measured 
with current technology. 

Model and methods .—We describe a 2DEC in a SC 
by means of a single-orbital tight-binding Hamiltonian 
of the form 

^ = E (4<^cj>+h-c.) ■ (1) 

(7j),cr 

This describes electrons with spin a =t, i hopping be¬ 


tween the nearest-neighbor sites (i, j) of a SC with Haus¬ 
dorff dimension dn = log^AT, as explained in Fig. 1. 
With reference to Fig. 1, each discretized SC sample is 
characterized by two integers, {n, m}. In this work n = 2, 
while m varies from m = 2 to m = 8. The Hamiltonian 
(1) is particle-hole symmetric and the spectrum of eigen¬ 
values extends from —4t to 4t for a bandwidth equal to 
8t. 

The hopping parameter t is used below as unit of en¬ 
ergy. Nanopattering a SC on the surface of semiconduc¬ 
tor hosting a high mobility 2DEC is expected [2, 7, 8] to 
yield t of the order of a few meV, while protocols based on 
STM manipulation [3] are expected to yield t ^ 100 meV. 
For the sake of simplicity, we ignore magnetic fields and 
electron-electron interactions, which are expected to lead 
to very interesting quantitative and qualitative effects 
that will be the subject of future works. 

The conductance and wavefunction calculations are 
performed by using KWANT [32]. In this toolkit, wave- 
function matching is implemented to compute the wave- 
functions in the scattering region and the scattering ma¬ 
trix Snm foi* an incoming propagating mode n and an 
outgoing mode m. The conductance between the left 
lead L and the right lead R is given by the Landauer 
formula 

G = x E (2) 

nGL,mGR 

The density of states (DOS) n{E) and optical con¬ 
ductivity cr(co’) are calculated by using the tight-binding 
propagation method (TBPM) [33, 34]. Since TBPM does 
not involve diagonalization of matrices, both CPU time 
and memory cost grown linearly with the sample size, 
allowing calculations with up to 10^^ sites. 

In TBPM, the DOS is obtained by the following 
Fourier transform [33-35] 

1 7^ 

t^{E) = ^ J exp {iEr/h) (<y?|<y?(r)) dr , (3) 

where the initial state \ip) = is a superposition 

of localized orbitals \i) with complex random coefficients 
Oi and the wave propagation = exp (—i'Hr/fi) \ip) 

is performed numerically by using the Chebyshev poly¬ 
nomial algorithm. Our results for the DOS of a {2,3} 
SC are reported in Fig. lb). We have checked that the 
function n{E) rapidly converges for a fixed value of n 
and increasing m: for example, results for the {2,7} and 
{2,8} SCs (not shown) are nearly indistinguishable over 
the entire bandwidth. 

Quantum transport and fraetal eonduetanee 
fluetuations. —We have calculated the energy de¬ 
pendence of the two-terminal conductance G{E) of 
the tight-binding model (1), for two different lead 
configurations, as in Figs. Ic) and d). 
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A summary of our main results for a {2,3} SC is re¬ 
ported in Fig. 2. In panel a) we clearly see that the two- 
terminal conductance G{E) is equal to /h for E = 0^ 
where a conductive extended state is present [17]. This is 
because for the two leads positioned as in Fig. Ic), elec¬ 
trons of a given spin injected on the left side of the sam¬ 
ple can reach the right side by following two equivalent 
paths, Fig. 2c), each carrying a conductance quantum, 
without being backscattered by the inner holes of the 
SC. On the other hand, as we can see from panel b), the 
sample can be insulating (i.e. G = 0), at the same energy 
when probed with the leads positioned as in Fig. Id). 

Data in Fig. 2 display also wild fluctuations. Such 
conductance fluctuations (CFs) can be quantitied by us¬ 
ing a box-counting (BC) algorithm [36]. This counts the 
number N of squares of size S, which is necessary to 
continuously cover the graph of G{E) (in units of e‘^/h) 
rescaled to a unit square. In general, points in the plane 
(log A^, — log (5) are expected to fall in three distinct re¬ 
gions. For large values of 6, the squares are too large to 
distinguish the features of the graph and N grows slowly 
as 6 decreases. For very small values of 6, the squares 
are so small that they resolve the single points in the 
set of data belonging to the CF graph: in this case N is 
expected to saturate to the number of points in the 
energy mesh where G{E) is evaluated. Finally, there is 
an intermediate region (usually called “scaling region”) 
where scaling is linear in the log-log plane, i.e. where 
N ^ 6~^. The slope d in the scaling region is the BC 
estimate of the Hausdorff dimension of the CFs. 

In Fig. 3, we show the results of the BC algorithm for 
the CFs of two SC samples with (slightly) different Haus¬ 
dorff dimensions, obtained by changing Af and C in the 
iterative geometrical construction illustrated in Fig. la). 
Results do not depend on the spatial arrangement of the 
leads. The analyzed CFs clearly show fractal behavior 
over a scaling region of more than two orders of magni¬ 
tude. Fractal CFs emerge [24] in the problem of phase- 
coherent ballistic transport and in the presence to strong 
quantum localization [36]. Indeed, in the SC transport 
problem, almost all wavefunctions are localized [28] [see 
e.g. Fig. 2d)] and do not contribute to transport (even 
in the absence of elastic disorder) due to scattering of 
electrons against the inner holes. 

Finally, we note that the value of d extracted from 
our numerical analysis is in excellent agreement with the 
Hausdorff dimension dn of the corresponding SC sam¬ 
ple [37]. It is remarkable that the analysis of CFs carries 
information on the sample geometry, down to very small 
length scales. 

Fig. 4 shows that the conductance of the 2DEG in a 
SC is robust with respect to both localized [panel a)] 
and smooth [panel b)] elastic disorder. The impact of 
a single-site vacancy—located on the highly-conductive 
path through the bulk of the SC displayed in Fig. 2c) —is 
shown in Fig. 4a). We see that, despite such a strong. 
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FIG. 2. (Color online) Energy dependence of the conduc¬ 
tance G{E) (in units of e^ jh) of a {2,3} Sierpinski carpet. 
Data in panel a) [panel b)] refer to the sample with leads 
positioned as in Fig. Ic) [Fig. Id)]. Panels c) and d) show 
the square modulus |'0(ri)|^ of the electron wave function 
corresponding to E = 0, with leads positioned as in pan¬ 
els a) and b), respectively. The colorscale varies from white 
= 0) to red (maximum). Black squares denote holes 
in the Sierpinski carpet. The wave function in panel c) rep¬ 
resents the state responsible for G(0) = 4e^//i in panel a). 
Note the extreme geometrical simplicity of this state, which 
connects the two leads while avoiding the holes in the Sier¬ 
pinski carpet. On the contrary, the wave function in panel d) 
represents an insulating state, where the electrons are totally 
back-scattered from the holes in the Sierpinski carpet. 

localized disorder source, G{E) still reaches its maximum 
value G{E) = ^e^ jh at E = 0. The impact of a smooth 
disorder potential, 

r = , (4) 

i 

which varies on an energy scale equal to 20% of the hop¬ 
ping amplitude, is shown in Fig. 4b). The spatial varia¬ 
tions of Sjii{ri)/t are reported in Fig. 4c), while the fate 
of the highly-conductive bulk path shown in Fig. 2c) is 
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FIG. 3. Box-counting analysis of the conductance fluctu¬ 
ations in the Sierpinski carpet. Data in panel a) have been 
obtained for a Sierpinski carpet with A/" = 8 and £ = 3, as in 
Fig. la). Data in panel b) have been obtained for a Sierpinski 
carpet with A/" = 12 and £ = 4. In each panel, + and x 
symbols represent the results of the box-counting algorithm 
applied to samples with the leads positioned as in Fig. Ic) 
and d), respectively. The horizontal dashed lines represent 
the constant value N = Ns, with Ns = 6000. In each panel, 
the slope d of the solid lines has been set equal to the Haus- 
dorff dimension dn of the corresponding Sierpinski carpet, 
i.e. d = 1.89 in panel a) and d = 1.79 in panel b). 


illustrated in Fig. 4d). 

Optical conductivity .—We have calculated the optical 
conductivity from the Kubo formula for the current re¬ 
sponse function to a uniform time-dependent vector po¬ 
tential. In the TBPM, this is calculated (omitting the 
Drude contribution at a; = 0) from [33, 34]: 


2 ^-/3huj 

Re[cr(u;)l = — lim - 

S e^0+ UJ 


-I 


sin (Hujt) 


x2 lm{(fi\f{H)J{T)[l - f{n)]J\(fi)dT ,(5) 

where S is the sample area, P = l/{k^T) is the inverse 
temperature, f{l-L) = 1/+ l] is the Fermi- 
Dirac distribution operator, and J(r) = ^lUrjhj^-lUrjh 
is the current operator in the Heisenberg picture of the 
time evolution. 

Fig. 5a) shows that the dependence of the absorption 
spectrum Re[cr(c(;)] (in units of /h) on the photon en¬ 
ergy huo (in units of t). This quantity is evaluated at 
a finite but small temperature /cpT t (specifically, 
/cpT = 0.0034t) and for two SC samples, i.e. {2, 3} and 
{2,8}. By looking at the difference between Re[cr(ci;)] 
at m = 7 (not shown) and m = S (and fixed temper¬ 
ature), we conclude that the data shown for m = 8 
represent the m ^ oo result. We note that the opti¬ 
cal absorption spectrum features three distinct behav¬ 
iors. For photon energies hw < 0.2t, absoption is very 
large, with Re[< 7 ( 0 ;)] ^ /h, peaking at « jh. The 

peak value is quite sensitive to the value of n in a {n, m} 
SC: for example, this peak reaches ~ AOe^/h in a {4,8} 
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FIG. 4. (Golor online) Quantum transport in a disordered 
Sierpinski carpet, with leads positioned as in Fig. Ic). Pan¬ 
els a) and b) show the energy dependence of the conduc¬ 
tance G(E) (in units of /h) in the presence of localized and 
smooth elastic disorder, respectively. Data shown in panel a) 
have been obtained by creating a single-site vacancy along the 
path of the conductive state shown in Fig. 2c), on the site of 
the Sierpinski carpet with spatial coordinates Vi = (10,18). 
Data shown in panel b) have been obtained by adding to (1) 
a smooth disorder potential Sja{ri) described by Hamiltonian 
(4). The potential is shown in c), where the colorscale varies 
from —O.lt (blue) to +0.1t (red). The conductive eigenstate 
at = —0.063t [indicated by a red arrow in b)] is shown in 
d). 


SC. For photon energies 0.2t ^ huj < 0.5t, Re[cr(cc;)] 
features a plateau at a value which is intermediate be¬ 
tween jh and lOe^/F. Finally, an exponential tail 
Re(cr) oc with a ^ 1.5 kicks in for huo > 0.5t. 

For a choice of the hopping amplitude t = 100 meV 
(which is the expected value for a SC realized by STM 
manipulation [3]), the absorption of the SC vanishes in 
the visible range of the electromagnetic spectrum and it 
is of the order of /h in the infrared. For photon fre¬ 
quencies spanning the range 300 GHz </< 5 THz the 
real part of the optical conductivity is very large, much 
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FIG. 5. Panel a) The real part of the optical conductivity 
a{uj) (in units of e^/h) of a 2DEG in a Sierpinki carpet is 
plotted as a function of photon energy huj (in units of t ), for a 
fixed value of temperature /cbT = 0.0034t (this choice implies 
T = 4 K for t = 100 meV). Data in this panel refer to two 
different values of m, m = 3 (thin dashed line) and m — S 
(thick solid line). Panel b) Re[cr(ct;)] (in units of /h) for 
m — t — 100 meV, and T = 4 K. 

larger than jh^ and displays a sequence of molecular 
resonances (whose smoothness is controlled by tempera¬ 
ture) . 

Summary .—In this work we have carried out exten¬ 
sive computer simulations of the conductance and op¬ 
tical absorption spectrum of a two-dimensional electron 
gas roaming on a prototypical plane fractal, the Sierpin- 
ski carpet. We have demonstrated that the conductance 
is sensitive to the spatial location of the leads and that 
it displays fractal fluctuations whose dimension is com¬ 
patible with the Hausdorff dimension of the sample. Due 
to the presence of a massive number of localized eigen¬ 
states, the optical absorption spectrum turns out to dis¬ 
play sharp molecular peaks at low photon energies, which 
may pave the way for applications in sensing and plas- 
monics. 
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